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ON FILTERED POLYNOMIAL APPROXIMATION ON THE 

SPHERE 

HEPING WANG AND IAN H. SLOAN 


Abstract. This paper considers filtered polynomial approximations on the 
unit sphere C obtained by truncating smoothly the Fourier series of 

an integrable function / with the help of a “filter” h, which is a real-valued 
continuous function on [0, oo) such that h(t) = 1 for t G [0,1] and h(t) = 0 for 
t > 2. The resulting “filtered polynomial approximation” (a spherical polyno¬ 
mial of degree 2L — 1) is then made fully discrete by approximating the inner 
product integrals by an A-point cubature rule of suitably high polynomial de¬ 
gree of precision, giving an approximation called “filtered hyperinterpolation”. 
In this paper we require that the filter h and all its derivatives up to 
are absolutely continuous, while its right and left derivatives of order 
exist everywhere and are of bounded variation. Under this assumption we 
show that for a function / in the Sobolev space VUp(S^), 1 < p < cxd, both 
approximations are of the optimal order in the first case for s > 0 and in 

the second fully discrete case for s > d/p. 


1. Introduction 

This paper is concerned with constructive polynomial approximation on the d- 
dimensional unit sphere C Our ultimate goal is an approximation scheme 

that uses only function values at selected points (sometimes called “standard in¬ 
formation”), but as an intermediate step we will consider also an approximation 
scheme that allows inner product integrals (“linear information”). The intermediate 
approximation scheme (defined precisely in Subsection 11.11 below) starts with the 
Fourier (or Laplace) series for an integrable function /, then truncates it smoothly 
with a suitable “filter” h. The final approximation is then obtained by approximat¬ 
ing the inner product integrals by a cubature rule of suitably high polynomial de¬ 
gree. Following [9], we shall call the resulting scheme “filtered hyperinterpolation”. 
We shall call the intermediate approximation “filtered polynomial approximation”. 

Remarkably, the convergence rate for both of the filtered approximation schemes 
is optimal (as explained in Remark 1.2 below or m), in a variety of Sobolev 
space settings. Precisely, under suitable conditions on the filter h and the cubature 
rule, and for s > d/p and 1 < p < oo, the Lp error of the fully discrete filtered 
hyperinterpolation approximation Vl.n (a spherical polynomial of degree 2L — 1 
defined in p.6p below) satisfies, for a function / in the Sobolev space Wp{S‘^) (see 
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Section [2] for definitions of the Sobolev spaces Wp{S‘^)), 

(1-1) \\f-VLMf)\\p<cL-’^\\f\\w^, 

where c depends on p, d, s and the filter h but not on / or L or the cubature 
rule. The condition s > d/p is the well known necessary and sufficient condi¬ 
tion for Wp{S‘^) to be continuously embedded in C(S'^), and is unavoidable in any 
approximation scheme that employs function values. The intermediate filtered ap¬ 
proximation scheme has the same rate of convergence, but no requirement on the 
smoothness order of the integrand / beyond s > 0. 

The optimal convergence result (HH) was proved by Mhaskar [5] in the context 
of neural nets for filters of C°° smoothness, and by Sloan and Womersley [5] for the 
case of p = oo and hlters of smoothness In the present paper we prove the 

result (lED for all values of p G [1,cxd], and for filters of smoothness = ^BV, 
where W^BV, r G N, denotes the set of all absolutely continuous functions h 
whose derivatives of all orders up to r — 1 are absolutely continuous and whose 

frl fr) 

right and left r-order derivatives h/_ \h_ oih exist everywhere and are of bounded 
variation, while [aj is the largest integer not exceeding a G R. For example for the 
important special case d = 2 it is sufficient that h G W^BV^ which certainly holds 
if h G For general d >2 we note that J C J BV and d -I- 1 > 

if d > 2, so we essentially improve upon the previous results. Moreover, we prove 
by counterexample that the present condition on the smoothness of the filter is not 
generally achieved for filters of smoothness merely ^ BV. 

The paper [8] constructed piecewise polynomial filters of smoothness merely 
C‘^~^ which generalise the de la Vallee-Poussin construction for d = 1, and studied 
their properties for the intermediate semi-discrete approximation and p = oo. Since 
these filters lie in W‘^BV, and since d > for all d > 2, these filters too are 

covered by the present theory. In the special case d = 2, even the de la Vallee- 
Poussin filter itself (which is simply a piecewise-linear (7° filter) is covered by the 
present theory. 

l.I. The approximation schemes. The Fourier series for / G is 

oo z{d,e) 

1=0 k=0 

Here (•, •) denotes the L 2 inner product 

(1-2) {f,g) ■■= f fgda, 

Js-^ 

where cr denotes the normalized rotationally invariant measure on S"^, Z{d,£) is the 
dimension of the space of homogeneous harmonic polynomials on R'^^^ of degree £, 
and {Yi^i.)e=o,i,...;k=i,...,z(d,e) is a complete set of orthonormal spherical harmonics 
on 

In the filtered polynomial approximation scheme the terms in the Fourier series 
are to be modified by multiplication by h{£/L), where h G C(R+) (a “filter”) 
satishes 

h{x) = 


1 if a; G [0,1], 

0 if a; > 2. 
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Thus the first stage of approximation is 

oo z(d,l) 

(1.3) / VlU) = Viif) := E E 

i=0 k=0 

Note that because of the properties of the filter we have VL{f) G well as 

Vkif) = / if / G n^, where 11^ is the set of spherical polynomials on of degree 
< L, which is just the set of polynomials of degree < L on restricted to S^. 
These are properties we want to preserve in the second stage of approximation. 

We obtain the second and final stage of the approximation by discretizing the 
inner product (11.21) by an iV-point cubature rule, 

N 

(1.4) / g{x)da{x) fv Qn{ 9) ■= '^Wig{yi), 

i=i 

where the points yi G and weights Wi > 0 are such that the cubature approxi¬ 
mation becomes exact for g G that is such that 

(1.5) [ P{x)d(7{x) = Qn{P) for P G n 3 ^_p 

(The 3L — 1 in this condition can be reduced to [(2-|-(5)L] — 1 for arbitrary but fixed 
<5 > 0, but for simplicity we shall stick with 3L — 1.) Thus the final approximation 
is defined by 

OO Z{d,l) 

(1.6) VlAJ) = V^E(/) := E E Hj:){f,Yi,^)Q.Yi,k, 

1=0 k=0 

where (•, •)q^ is the discrete version of the inner product (EH), defined by 

{f,9)QN ■= QN{f9)- 

Subsequently, the notation A{N) x B{N) means that A{N) <C B{N) and 
A{N) ^ B{N), and A{N) ^ B{N) means that there exists a positive constant c 
independent of N such that A{N) < cB{N), and A{N) ^ B{N) means B{N) <C 
A{N). We write also a+ := max{a,0}. 

The main result of the paper can be formulated as follows. 

Theorem 1.1. Let h he a BV filter, and let Vl and Vl,n be as in dOl) 

and (fTel) with dig respeetively. Then for all f G Wp{Ei‘^), 1 < p < oo, s > 0, the 
inequality 

(1-7) \\f-VL{f)\\p^L-^\\f\\wi 

holds, and if s > d/p then (11.11) also holds, that is 

(1-8) ll/-VL,iv(/)||p«T-1l/lk- 

Remark 1.2. For iV G N, the sampling numbers (for optimal recovery) of the 
Sobolev classes BWp{S'^) (the unit ball of the space Wp{S'^)) in Lp{S‘^) are defined 
by 

gM{BW;{^<^),Lp{^‘^)):= inf sup \\f - ^{f{^i),..., f{^N))h, 
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where the infimum is taken over all N points .. ,^n in and all mappings (p 
fromR^ to Lp{S‘^). It follows from [13] that 

(1.9) gN{BW;iS‘^),Lp(S'^))-N-^/‘^. 

To express the error for filtered hyperinterpolation in terms of N rather than L, 
we use the known facts that the minimal number of points N for a cubature rule 
with polynomial degree of precision L is N (see [1], Theorem 7.1), and that 

cubature rules with the number of points N of this order do exist (see M)- For such 
optimal cubature rules it follows from (HD that 

sup \\f-VL,N{f)\\p X X gN{BW;{S‘‘),Lp{n‘^)), 

which implies that the filtered hyperinterpolation operators Vl,n OLre asymptotically 
optimal algorithms in the sense of optimal recovery. 

The paper is organized as follows. Section 2 contains some basic notation and 
properties concerning harmonic analysis on the operators Vl and Vl,n- In 
Section [3] we prove Theorem 1.1. The proof is modelled on one in [12] for the 
hyperinterpolation approximation, where hyperinterpolation is the approximation 
obtained by replacing the filter h in p.6ll by the characteristic function of [0,1] (see 
0)- 


2. Preliminaries 


Let = {a; S : |a;| = 1} (d > 2) be the unit sphere of (with 

I a; I denoting the Euclidean norm in endowed with the rotationally invariant 

measure a normalized by da = 1. For 1 < p < oo, denote by Lp = Lp{S‘^) the 
collection of real measurable functions / on with finite norm 

ll/llp = \f{x)\P da{x)^ " < +CX), 1 <p < cx). 


For p = oo, the essential supremum is understood instead of the integral. 

We denote by TLf the space of all spherical harmonics of degree I on It is 
well known that the dimension of TLf is 


Z{d,£) := dimH 



1 , 

(2f+d-l) {e+d-2y. 
(d-1)! V. 


if £ = 0, 
if £ = 1 , 2 ,..., 


and 

Z{d,£)^l‘^-\ £ = 1,2,.... 

Let 

{Ypk = Y,% \k = l,...,Zid,£)} 
be a fixed orthonormal basis for TLf. Then 


{Ypk \ k = l,...,Z{d,£), £ = 0,1,2,...} 

is an orthonormal basis for the Hilbert space L 2 {S‘^). Thus any / G L 2 {S‘^) can be 
expressed by its Fourier (or Laplace) series: 


cxD OO Z{d,l) 

f = Y,He{f)=Y^ ^ {f,Yp,)Ypk. 

i=o e=o k=i 
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where He{f) = £ = 0,1,, denote the orthogonal projections of / onto 

T-Lf, and 

{f,Ye,k)=[ fix)Ye^kix)da{x) 

are the Fourier coefficients of /. Let Cg{t) be the usual ultraspherical (Gegenbeuer) 
polynomial of order A normalized by C'^(l) = r( 2 i^)r^i+i) generated by 

1 °° 

+ <.<!)■ 

(See ([in], P- 81). The operator He can be expressed as follows: 

z(d,e) 

Hef{x)= ^ {f,Ye^k)Ye^k{x) = 

k=l 

where 

z{d,e) 

Ke{x,y)= Yi^k{x)Ye,k{y) = ^^Ce{x-y), a:, ?/e A = —d > 2, 

fc=i 

which is a “zonal” kernel, that is one depending only on x • j/. 

It is well known that the spaces Hf, £ = 0,1,2,..., are just the eigenspaces 
corresponding to the eigenvalues —£{£ + d— 1) of the Laplace-Beltrami operator A 
on the sphere. Given s > 0, we define the s-th order Laplace-Beltrami operator 
(—A)'* on in a distributional sense by 

Hoii-Anf)) = 0, Heii-Anf)) = {£{£ + d- l^Heif), £=1,2,..., 

where / is a distribution on For f & Lp, if there exists a function g £ Lp such 
that He{g) = {£{£ + d — l)y/'^He{f) for all £ e N, then we call g the s-th order 
derivative of / in the sense of Lp and write g = = (—A)®/^/. For s > 0 and 

1 < p < oo, the Sobolev space Wp = Wp{S^) is defined by 

W;:={f£Lp : ll/lk. :=||/||p+||/(^)||p<(X)}, 

while the Sobolev class BWp is defined to be the unit ball of Wp. For / £ Lp, 1 < 
p < oo, we define 

EL{f)p-.= M{\\f-P\\p:P£ni}, L£N. 

It follows from the Jackson inequality (see [Tl]) that 

Az.(/)p«L-ni/|k- 

Let h be a “filter”, i.e., a continuous function on [0, oo) such that h{t) = 1 if 
t £ [0,1] and h{t) = 0 if £ > 2. For L > 1 we define a zonal kernel by 

OO ^ 

$L(a:, y) = $L(a:, y) =Y Kj)Ke{x, y), x,y £ 

i=0 


f{y)Ki{x,y)da{y), 


( 2 . 1 ) 
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and the operators Vl by 

oo Z{d,£) oo j 

e=o k=i e=o 

(2.2) =[ fiy)^L{x,y)da{y), f G Lp , 1 <p < oo. 

It follows from the definition of Vl that Vl^/) G ^ 2 L-i Yl{P) = P if P S 
For a linear operator A on Lp{S‘^), 1 < p < oo, the operator norm ||A||(p_p) of A 
is defined by 

P||(p,p) :=sup{P/||p :/GPp(S"), ||/||p<l}. 

For simplicity we write ||A|| = ||^||(oo,oo)- Using a duality argument we obtain 

(2.3) ||Ul|| := IIVlII(oo.oo) = I|Ul||(i,i) = sup f \<t>L{x,y)\da{y) < oo, 

xOS’‘ JS'^ 

noting that the last integral is independent of x because of the zonal nature of $. 
It follows from the Riesz-Thorin theorem that for I < p < c», 

linii(p.p) < (|inii(i,i))'/'’(||yL||(oo.oo))'”'^'’ = IIVlH < cx). 

Since the operators Vl are linear operators on Lp{S‘^) satisfying Vl[P) = P for 
P G n^, by the Lebesgue theorem the error of approximation by Vl in the Lp 
norm can be estimated for s > 0 as 

11/ - Ul(/)|Ip < (1 + ||Ul||(p.p))Pl(/)p < (1 + ||Ul||)Pl(/)p 

(2.4) «(l + l|UL||)U-1|/|k«, /GIU;(S"). 

The last inequality in p.4ll shows that estimation of the operator norm ||Ul|| is 
extremely important in numerical computation. We now turn to this question. 

A function / on [a, b] is of bounded variation and denoted by / G BV[a, b], if 

n 

WJWbv ■■= sup'^\f{xk) - f{Xk-l)\ < oo, 

where the supremum is taken over all partitions T : a = xq < xi < ■ ■ • < Xn = b. 
As above, for r > 1 we denote by BV[a,b], r G N the set of all absolutely 
continuous functions h on [a,b] for which is absolutely continuous and 

and exist and are of bounded variation on [a, 5]. It is easy to verify that the 
function f{t) = (1 — t/2)“, a > 0 belongs to IUL“JPU[0,3] and does not belong 
to cU+iJ [0, 3]. The following lemma gives a condition on h for which the operator 
norms ||Vl|| are uniformly bounded, i.e., 

sup IIVlII < oo. 

L>1 

Lemma 2.1. Let h be a W'^^'^BV filter and let Vl be given as in (12.211 . Then 
we have 

(1) for any f G Lp(S‘'), 1 < p < oo, Ul(/) G Ul(P) = P for all 

P G ni; 

(2) the operator norms ||Vl|| are uniformly bounded; 
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(3) for f e Lp(S^), l<p<oo, 

I|Vl(/)|Ip«||/|Ip; 

(4) for any f G Lp(S‘'), 1 < p < cx), 

||/-VL(/)||p«i?L(/)p; 

(5) for any f G Wp(S‘^), I < p < oo, s > 0, 

||/-^L(/)||p«L- 1 |/|k- 

Proof. It suffices to prove (2). Given 6 > —1, the Cesaro (C, ^)-operators of the 
spherical harmonic expansion are defined by 

1 ^ 

<^Uf) = ^T.^i-kHkif), L = 0,1,2,..., 

^ k=0 

where Af. = fc = 0,1, 2,.... It is well known that the operator norms \\crj^\\ 

of the Cesaro operators (t£ are uniformly bounded if and only if 5 > (see (TU 
p. 55]). 

Let / be a real-valued function defined on R. For r = 0,1, 2,..., we define the 
difference operator A*" by 

A0 = /, Afix) = fix)-fix + l), A’-/(x) = A(A’'-V(a:)), 


where I denotes the identity operator. By induction, it is easy to see that 

k=0 ^ ' 

For a sequence the difference operator A’'afc is defined as A’'/(fc) with 

/(fc) = Ofe, fc G N. If h G W^'BV, then we easily verify that 

A^h{y) = {-lYL-^ [ h^Y\^±ffl±-A±ffL)du,---dur. 

For r = , the summation by parts formula applied to (12.21) gives 


2L 


(2.5) 


Jr Jr 

VUf) = A^^^h{-)AlaUf) = E ^^^^h{-)Alal{f). 


/c=0 fc=0 

(See [21 pp. 408-409].) For a W^BV filter h we have 

21/ , 2L ^ _i_ 1 Z- 

^ \A^+^h{j)\ = ^ lA’-h(^) - A’'h(-)| 






^L' 


2L 


= r-''Y.\ 


k=0 


< L' 




[ 0 . 1 ] 

2L 

AI A( 


(r) ,k + 1 + Ui + ■ ■ ■ + Ur . ,(r),k + Ui + ---+Ur., 

hyi - - -) - hyi - - -) I dui 


(r) / fc + 1 + Wl + ■ ■ ■ + Ur , Ur , 

‘ - j -) - k\\ - - -) 


= 

k^(/)lloc« 


We note that for r = and fc = 0,1,..., 2L, 


and |^^|<(fc + l)’-<(2L)L 


■ dUr 


dui ■ ■ ■ du. 










HEPING WANG AND IAN H. SLOAN 


It then follows from (12.511 that 

2L , 

rL(/)iioo < E Ai « L-’'(2L)ni/iioo«ii/iloo, 

k=0 

which completes the proof of Lemma 2.1. □ 


The following lemma gives an example of a filter h of smoothness merely ^ BV 

such that the corresponding operator norms || V^|| are not uniformly bounded. This 
implies that the optimal order of convergence is not generally achieved for filters 
of smoothness merely ^BV. So the condition on the smoothness ^ BV 

of the filter h for which supj;, II < oo is tight. On the other hand, we note that 
this does not preclude the possible existence of particular filters of smoothness 
merely ^BV for which sup^ ||b^L II is finite. 


Lemma 2.2. There exists a filter h of smoothness merely hhi 2 \bV such that 


sup||V£‘|| = 00 . 


Proof. Let h^ft) := (1 — t/2)^ ^ i^g g (^00 gj^gj. gg [0,oo) such that 

hi[t) if t S [0,1] and hi{t) = 0 if t > 3/2. We set 

h{t) := ho{t) - ho{t)hi{t) + hi{t) = ho{t) (1 - hi{t)) + hi{t). 

Obviously, h satisfies h{t) = hi{f) = 1 if t G [0,1] and h{f) = ho{t) if t > 3/2, 
thus h is a filter, but it is a filter of smoothness merely ^BV. We note that 

h 2 {t) := —ho{t)hi{t) + hi{t) belongs to (7°°, so the corresponding operator norms 
||y^^|| are uniformly bounded. On the other hand, when h{t) = (1 — t/2)^, the 
corresponding operators are the Riesz operators R^l- By the equivalence (due 
to Marcel Riesz) of the Riesz and the Cesaro summability methods of order d > 0 
(see 0). we deduce that the Riesz operators are uniformly bounded if and only 
if d > It follows that the operator norms ||V^‘’|| are not uniformly bounded, 
so the operator norms WVjfW = ||V^° + V^^|| are also not uniformly bounded. □ 


In the sequel, we always assume that h is a, ^ ^ BV filter. It follows from 
(12.31) and Lemma 2.1 that 

(2.6) ||Vl|| = sup ||$l(x, Olli < 1- 

For / G IFp (S'^), 1 < p < 00 , s > d/p, we define 

ro{f) = Vi{f), Tr{f) = V 2 r{f)-V 2 r-.{f) for r > 1. 

00 

Then Tr{f) G Il^r+i, and ^ Tr{f) converges to / in the uniform norm. Further- 

r—0 

more, by (12.41) and (12.6p we have 

(2.7) ||t,(/)||p < 11/ - V2rif)\\p + 11/ - R2.-i(/)||p « 2—ll/lliy^. 
for s > 0. 
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For L e N, we assume that QN{g) ■= is a positive cubature rule 

on which is exact for g G i.e., the points yi G S'^, weights Wi > 0, satisfy, 

for all g G 

N 

(2.8) / g{x)da{x) = Qwig) ■= '^Wigijji). 

i=i 

The filtered hyperinterpolation operator Vl,n is dehned by 

oo Z{d,£) 

VLMf) = vtN{f) :=E E 

£=0 k =0 

where (•,-)qjv is the discrete version of the inner product, defined by {f,g)Qf^ '■= 
QNifg)- According to (I2.ip . we have 

N 

(2.9) Vl,nU){x) = {f,^L{-,x))Qr, 

f=l 

Note that for any P G 11^ we have P{-)^l{x, •) G for arbitrary fixed x G 

Thus by (12.8p we get 

(2.10) P{x) = Vl{P){x) = {P.^l[-.x)) = {P,<^l{:x))q^ = Vl,n{P)- 
We are now ready to prove Theorem 1.1. 

3. Proof of Theorem 1.1 

In order to prove Theorem 1.1, we need the following lemma. 

Lemma 3.1. [2j Theorem 2.1] Suppose that T is a finite set ofS‘^, and : w G T} 
is a set of positive numbers satisfying 

E « / \Pix)rdaix), P G nt 

<.;er 

for some 0 < po < oo and some positive integer n. // 0 < pi < oo, m > n, and 
P G then 

(3.1) E « i-r [ \p{y)rdaiy). 

Proof of Theorem 1.1. The inequality (E3 follows directly from Lemma 2.1 (5), 
so it suffices to prove m- First we consider the case p = oo. Equation (12.101) 
combined with the Lebesgue theorem implies that for / G VF^^PS^^) and s > 0 

(3.2) 11/ - 14,w(/)||oo < (1 + \\VL,N\\)EL{f)oo « (1 + \\VL,N\\)L~^\\f\\w^- 
Using the standard argument, we obtain from (1211) 

N 

\\Vl,n\\ = sup E^jl^i(a;.2/i)l- 

By ()2.8I) . we know that the set {j/i,..., j/jv} satisfies the condition of Lemma [3.11 
with pq = 2. It then follows from Lemma l3 .1 1 with pi = 1 that 

II Vl,ATII < sup 



\‘^L{x,y)\da{y) <C 1, 
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where in the last step we used (|2.6I1 . Then (|1.8I1 for the case p = oo follows directly 
from (j3.2|l . 

Now we consider the case 1 < p < oo. We noted already that for / G s > 

OO 

d/p, the series ^ Tr{f){x) converges to / uniformly in Given L > 0, let m be 

r—0 

the integer satisfying 2™ < L < 2"*+^. Since Tr{f) € we get from (12.101) that 

for / G with s > p/d, 


OO 

f{x) - VL,N{f){x) = {Tr{f){x) - VL,N{Tr{f)){x)) . 

r—m 


It follows that 


11 / - Vlw(/)IIp < E - ^Giv(T-.(/))||p 

r^m 

(3.3) < E {\\rrfh+\\VL,N{rr{f))h). 

r^m 


We claim that for 1 < p < oo and r > m. 


(3.4) 


l|VLw(r.(/))||p« 



P- 


In fact, for p = 1 and r > m, by (j2.6D we obtain 


N 

||nw(r.(/))||i = |lE^i^-(/)(%)‘fL(-,%)lli 

N 

i-i 

N 

«(^)'ik.(/)iii, 


where in last step we used again Lemma 3.1 with pi = 1. 
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For 1 < p < oo and r > m, by the Holder inequality and (|2.6|1 we have 



— 1 /* 

< (max||$z,(x, Olli)^ \Trif){yj)\P \^L{x,yj)\^da{x) 


N N 

|T.(/)(y,)ni<i>L(-,%)||i \rr{f){yjW 

3=1 i=i 



where we twice used Lemma 3.1, with pi = 1 and pi = p. Thus (|3.4p is proved. 

It then follows from (j3.3p . (13.411 . and (j2.7p that for / € IT*1 < p < oo, and 
s > d/p, 

ii/-u.^(/)iip« E + 

r>m 

^ 2-md/p 

r>?7i 

which proves (11.81) for 1 < p < oo. The proof of Theorem 1.1 is now complete. 
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